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In  1973  Lea  and  Boncal  Cl]  converted  the  Oceanics  Inc-  6 
DOF  Loads  and  Motions  digital  computer  simulation  program 
for  the  100-B  Surface  Effect  Ship  to  reflect  the  physical 
characteristics  of  the  X8-3  craft.  Since  tnat  tine  there 
has  been  a continuous  effort  at  the  N PS  to  update  and 
improve  on  this  simulation  [2, 3,4,5]. 

One  of  the  major  drawbacks  of  that  program  which 
precludes  regular  use  of  it  for  sea  state  operations 
simulation  is  t.ne  excessive  computational  and  turn-around 
times  required  [6]  . As  it  is  presently  desired  to  use  the  6 
DOF  simulation  for  design  and  evaluation  studies  of  the 
proposed  3K  Ton  5SS  control  System  in  varying  sea  state 
conditions,  it  is  important  to  significantly  reduce  the 
simulation's  computational  time. 

One  area  where  the  possibility  exists  to  achieve  a major 
reduction  in  computational  time  is  tie  manner  in  which 
sidewall  forces  and  moments  are  calculated.  The  present 
form  of  the  6 DOF  program  uses  a large  number  of  discrete 
sections  (28  per  sidewall)  with  an  off-line  table  look-up 
scheme  to  determine  the  forces  generated  in  each  discrete 
section.  The  simple  model  derived  by  McIntyre  C7]  uses  two 
sections  per  sidewall  with  very  simplified  geometry  that 
permits  on-line  computation  of  the  forces  and  monents.  The 
result  of  eventual  incorporation  of  tais  sidewalL  geometry 
and  forces  calculation  into  the  5 DOF  program  should  be  a 
significant  reduction  of  computation  and  turn-around  time. 

Another  area  of  great  interest  in  SES  Simulation  is 
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ta  a 


gaining  insight  into  tae  dynamic  behavior  of 


craft ' s 


vertical  motion  of  pitch  and  heave  due  to  pressure  induced 
lift  forces  generated  in  the  plenum  and  seals.  As  stated  in 
ref  [3]  as  recently  as  August,  1 973, 

•'The  overall  precision  of  the  analysis  is  satisfactory 
for  predicting  stability  and  maneuverability.  However, 
several  areas  present  vecy  complex  or  difficult  physical 
pnenomena  to  analyze  and  further  effort  should  be 
invested  to  improve  understanding  of  tnese. ' ' 

The  first  area  for  continuing  investigation  reconaended  was; 

' 1 A oetter  determination  of  pitch-heave  dynamic 
cnaracteristics.  • ' 

An  improved  knowledge  of  tne  dynamic  characteristics 
brougnt  about  by  the  pressure  induced  forces  and  moments 
could  lead  to  a reduction  of  the  computation  time 
required  for  craft  simulation.  These  pressure  effects 
introduce  a stiffening  in  the  integration  of  sea  state 
induced  oscillations  and  a better  understanding  of 
processes  could  lead  to  the  removal  of  this  stiffness 
and  thus  reduced  TPtI  time. 


In  this  work  both  of  these  areas  are  considered.  First 
of  all  the  various  craft  physical  characteristics  are 
modeled,  resulting  in  Changes  and  improvements  to  the 
initial  efforts  by  McIntyre  C7]  . Tnese  equations  are 
linearized  about  a steady-state  operating  point  using 
the  Taylor  Series  expansion.  The  linearized  equations 
are  cast  into  a State  Space  format  and  simulated  on  the 
digital  computer.  loraparison  is  made  between  the 
time-domain  results  of  the  linear  model  and  tne  6 DOF 
model.  The  State  Space  representation  is  used  to 
generate  a signal  flow  graph.  Mason’s  Gain  Rule  is 


applied  and  the  S-iomain  roots  of  the  craft’s 
characteristic  equations  are  located  and  it  Is  shown 
that  soma  simplification  and  redaction  of  system  order 
is  possible.  The  resultant  S-iomain  polynomials  are 
used  to  generate  Bode  Plots  to  graphically  display  the 
craft's  simplified  frequency  response. 

Conclusions  are  drawn  and  recommendations  are  made 
concerning  the  validty  of  the  model  and  areas  for  future 
study. 

Several  procedural  notes  are  in  order.  Because  all 
equations  are  to  be  solved  digitally  using  the  FORTRAN 
computer  language  and  because  dijital  computer  text 
processing  facilities  were  used  to  generate  this 
document,  all  formulas  will  be  presented  in  FORTRAN  I 

format.  in  the  thesis  oody  reference  to  the  ' 1 nonlinear 
model' ' , ''six  degree-ot  -freedom  model''  or  ' '6  DOF'' 
model  ail  mean  the  Loads  and  notions  program  as  adapted 
to  the  XR-3  test  craft  at  the  Naval  Postgraduate  School 
[3]  . Reference  to  the  linear  model  refers  to  the 
linearized  model  of  the  same  craft  as  developed  herein. 
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A.  BACKGROUND 


The  original  development  wart  on  the  simplified 
two-deg ras -of  -freedom  model  was  done  oy  McIntyre  C7J  . 
Several  weaknesses  i.a  the  original  model  wera  noted  by 
McIntyre  in  that  work. 

The  major  portion  of  the  work  presented  here  was 
involved  in  close  examination  of  all  aspects  of  the 
original  model,  comparison  of  tne  modeling  techniques 
used  therin  with  other  methods  and  remodeling  of  several 
important  craft  pnysical  char  act ar i s t ics . 

One  major  area  of  invast igat ion  was  tne  lack  of  pitch 
damping  in  the  original  model,  which  resulted  in  the 
addition  of  an  added-aass  effect  to  the  pltca  moment 
equat ions . 

Additional  changes  involve  a revision  of  the  force 
equations  for  the  rear  seal,  major  changes  ia  the  way 
tne  plenum  chamoer  is  modeled  (also  added  to  tae  loads 
and  motion  o-dof  model  [23  ) . and  a more  rigorous 
development,  of  the  planing  force  a (nation  of  [7]  . 

There  are  several  areas  of  tneory  concerning  CAS  3ES 
modeling  where  empirical  data  is  lacking  due  to 
difficulties  with  physical  measursmen cs.  One  such  area 
inviting  closer  inspection  was  dynamic  variation  of  the 


effective  center  j£  pressure  Cor  the  plenum  gauge 
pressure  lift  effect.  \ possible  approach  is  suggests! 
herein  and  included  in  the  model. 


B.  ASSUSPTIONS 

The  following  assumptions  are  made  for  this  constant 
speed,  pitch-heave  model: 

1.  Planing  forces  are  concentrated  at  tha  centroid  of 
the  keel  longitudinal  cross  section. 

2.  Total  seal  force  is  Lumpel  into  a pressure 
differential  acting  on  the  seal  face  area  in  contact 
with  the  water  passing  beneath  the  seal. 

J.  Aerodynamic  forces  on  the  forward  seal  are 
disregarded. 

h.  Dynamic  variations  in  moment  lever  arms  are 
disregarded  for  all  moments  except  seals  and  the  moment 
due  to  plenum  guage  pressure. 

5.  Dynamic  variations  in  effective  plenum  roof  area 
will  be  included. 

a.  Any  lift  effects  lue  to  aerodynamic  forces  are 
Lumped  in  with  the  planing  forces  (for  reasons  to  be 
ieveloped  in  Ihapt  IV.)  The  moments  are  included  in 
plenum  gauge  pressure  moments  . 

7.  Constant  seal  leakage  area  is  assumed. 

3.  Since  this  aodeL  is  for  constant  surge  velocity. 
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drag  and  thrast  forces 
equilibrium  and  in  balance 


are  considered  to  be  in 
for  moments  about  the  Y axis. 


C.  COORDINATE  SYSTEM 


(see  fig.  1 ) 


1.  The  origin  of  tae  cartesian  coordinate  system  is 
located  at  the  calm-water  line,  longitudinally  and 
transversly  coincident  with  the  craft  center  of  gravity. 

2.  The  X axis  is  the  longitudinal  coordinate  defined 
as  positive  as  X increases  towards  the  bow  of  tae  craft. 


2.  The  Y 
positi vely 


axis  is 
towards 


the  transverse  coordinate 
the  starboard  side. 


increasing 


4.  The 
positively 


Z axis  is  the  vertical 
in  a downward  direction. 


axis. 


increasing 


5.  ZS  is 
gravity  to 


the  vertical 
the  keel. 


distance  from 


the  canter  of 


5.  Z is  the  vertical  distance  from  the  calm  water  line 
to  the  vertical  center  of  gravity. 

7.  THETA  is  the  pitch  angle  defined  as  positive  as  the 
bow  pitches  up. 


3.  Positive  pitch  moments  are  defined  as  moments  which 
tend  to  cause  a positive  increase  in  the  angle  THETA. 


D.  SIDEWALL  MODELING 


15 


As  it  was  desired  to  keep  ail  physical  dimensions 
consistant  with  the  six  decree  of  freedom  modal,  scale 
model  drawings  ware  developed  with  the  data  from  the 
sidewall  subroutine  of  that  model,  All  sidewall 
dimensions  usad  in  this  aodel  were  taken  from  those 
drawings. 

The  sidewalls  are  modeled  in  two  sections,  ona  forward 
of  tne  longitudinal  3.3.  and  the  other  aft.  Each 
sidewall  is  considered  to  be  of  constant  width  at  the 
water  line  and  the  keel,  with  width  increasing  in  the 
vertical  direction  according  to  an  average  deadrise 
angle. 

The  average  forward  sidewall  area  is  found  to  be: 

(see  fig.  2 and  3b  ) 

AV2RAG2  < IDTH  * LD  B A 3 1 / ( 2 * T AM  pH  1 ) ) ► wS  1 0 


(1 1-1) 

where  LDBARl  is  draft  at  tocwacd  centroid,  DR  1 = 
average  deadrise  angle  forward  and  M 5 1 0 is  the  average 
keel  width  of  the  forward  section.  Therefore  submerged 
volume  of  a forward  sidewall  is  given  oy: 


VOL  = LI* (LDBAR1/(2*TAN ( OR  1 ) ) ♦HSIO) *LDBAR1 


(II-2) 


where  LI  is  the  length  of  a forward  ouoyant  section. 
Following  the  same  procedures  the  after  sidewall  width 


lb 


— 


1 


is: 


(sae  fig.  2 aa d 3c  ) 


A V 2RAGE  WIDTH  * LD8A 52/ ( 2 *TAN (DR2U  * WS2Q 


(H-3) 

where  DR2  * average  deadrise  angle  aft 
LD3AR2  = draft  at  after  buoyant  centroid 
WS2Q  3 average  iceel  width  aft 

Submerged  volume  of  an  after  sidewall  is  then 
determined  to  be; 


VOL  3 L2*  (WS20  + LD3A2  2/(2*TAN(DR2) ) ) * LD3 A 32 

(II-<0 

where  L2  is  the  length  of  the  after  sidewall. 


S.  SEAL  MODELING 


A major  change  in  dclntyre's  model  was  made  by 
remodeling  the  stern  seal.  In  the  original  model  a large 
negative  moment  was  required  to  place  the  pitch 
aquations  in  steady  state  equibibriua.  Investigation  of 
seal  forces  as  presented  in  R]  revealed  a large 
discrepancy  between  the  after  seal  forces  generated  in 
Sclntyre's  model  and  calculations  based  on  the 
information  in  GO  . This  discrepancy  was  further 
supported  by  actual  seal  load  data  presented  oy  Layton 
in  Cldl  . All  seal  dimensions  used  in  this  model  were 
taiten  from  data  presented  in  [11]  . Seal  hinge  Locations 

were  taken  from  the  data  used  in  G3  . 
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Seal  lift  force  forward  is  generated  Dy  plenum  gauge 
pressure  acting  downwird  on  the  seal  face  area  in 
contact  with  the  water  passing  underneath  the  seal. 
This  downward  force  results  in  a lift  force  being 
generated  as  the  downward  movement  of  the  seal  displaces 
the  water  beneath  it. 

The  seal  is  considered  to  have  a straight  forward  face 
which  runs  from  the  hinge  point  at  the  plenum  roof  to 
the  oottom  of  the  'teel  at  an  angle  determined  to  be 
best-fit  from  the  data  in  ref.  DO  • The  wetted  surface 
of  the  seal  is  considered  to  lay  parailei  to  and  on  the 
water  at  the  point  where  the  seal  face  intercepts  the 
calm  water  line,  as  a result  forward  seal  lift  equations 
are: 

ASEAL1  * WIDTH*  XSEAL  1 

(II-5) 


XS2AL1  = (LD-L3*TAM (THBTA)  ) /SIN (31) 


( I 1-6) 

LD  - draft  at  the  2.3 

L.3  * distance  from  2.3  to  point  of  contact  for  forward 
seal  face. 

LD  - L3*TAM  (THETA)  = draft  at  forward  seal  face  contact 
point 

31  degrees  is  the  interior  angle  of  the  seal. 

The  after  seal  wetted  area  is  determined  in  a similar 
manner  to  be  : 

A3EA L2  = W IDT H* XSE A L 2 
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(II-7) 


XSSAL2  = ( LD*L4  *IAN (THETA) ) /SIN (32) 

(II-9) 

In  this  modal  forward  seal  pressure  is  considered 
equal  to  plenum  pressure.  As  a result  forward  seal  lift 
is  : 


H52ALF  = PB3AH* AS2AL1 


(11-4) 


Since  plenum  guage  pressure  tends  to  force  the  after 
seal  upward,  resulting  in  venting  of  the  plenum  air 
mass,  the  rear  seal  is  operated  at  a pressure  slightly 
above  that  of  the  plenum  to  reduce  the  venting  tendency 
of  the  after  seal.  As  a result  taa  after  seal  requires  a 
modification  in  modeling.  In  dclntyre's  work  the 
pressure  differential  (PSE AL-P5BAR)  only  was  considered 
in  the  lift  force.  However,  as  previously  stated,  there 
were  large  discrepancies  in  computed  loads  as  compared 
to  measured  data  from  [10]  . Therefore  in  this  model  , 
the  seal  woricing  pressure  is  considered  to  be  ? 3 3 A R plus 
differential  pressure.  As  a result  after  seaL  lift  is 
given  by: 

H3SALA  * (P3BAR  + PDIFF)  *A5EAL2 


(11-10) 


F.  PLENUM  MODELING 

la  the  wortc  preseated  Dy  McIntyre  the  planus  was 
considered  a box  with  the  oniy  dynaaic  dimensional 
variation  being  in  the  Z direction.  However  (as  can  be 
seen  in  fig.  1 ),  although  the  forward  edge  can  be 

assumed  to  be  reasonably  constant  in  tne  Y-2  plane,  the 
forward  face  of  the  rear  seal  slopes  rather  steeply. 
Therefore  the  plenum  Length  (in  the  !C  direction  varies 
rapidly  as  draft  varies.  As  a result  the  effective 
plenum  roof  area  consists  of  that  area  from  the  after 
edge  of  the  forward  seal  to  the  forward  hinge  point  of 
the  rear  seal,  plus  the  vertical  projection  of  the 
unwatted  forward  face  of  the  raar  seal.  As  an 

additional  consequence  it  can  be  seen  that  tne  effective 
plenum  center  of  pressure  varies  fore  and  aft  as  draft 
varies  . In  fact  as  draft  increases  the  Z.P.  moves 
forward. 

The  following  variable  names  will  oe  used  to  describe 
the  plenum  chamber  (see  fig.  u . ) 

L?3  - length  of  plenum  roof  from  the  rear  elg?  of  the 
forwari  seal  to  the  hinge  point  of  the  after  seal. 

L?'4L  3 length  of  plenum  at  the  heel  from  the  after  edge 

of  the  forward  seals  to  the  transom. 

LP  * length  of  plenum  at  the  waterline. 

3U3HGT  * vertical  iistince  from  the  plenum  roof  to  the 
bottom  of  tne  heel  at  tne  center  of  gravity. 

AB  * plenum  effective  roof  area. 

Using  the  previous  lefinitions,  tne  plenum  roof  length 
and  area  is: 
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L?  = LP2L-  (LPWL-LPR)  *(LD/3UBHGT) 


and  effective  plenum  roof  area  is 


(II-  1 la) 


AB  = aiDTH*LP 


(II-  1 1 b) 

with  the  result  that  lift  due  to  planum  gauge  pressure 
is: 


H ? 3 ES  = -A3*PBBAR 

(II-11C) 

As  a comment  on  tne  complexity  of  the  lift  forces  of 
this  model  note  that  for  constant  plenum  pressure 
increased  draft  results  in  decreased  plenum  lift. 

G.  CRAFT  DYNAMICS 


1.  Plenum  Air  Mass  and  Pressure  Equations 

For  the  purposes  of  this  modal  tne  raar  seal 
leakage  area  is  considered  constant  since  to  assume 
otherwise  would  lead  to  modeling  a very  difficult  and 
poorly  understood  dynamic  process.  In  addition  the 
plenum  pressure  VS  air  mass  process  is  assumed  to  be 
adiabatic  in  nature.  Fan  characteristics  are  tnose  used 
in  the  constant  air  leakage  heave  only  model  of  Gerba 
and  Thaler  in  Du3  . Air  flow  into  the  plenum  is  taken 
as  positive. 
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herefore  the  airflow  into  the  chamber  is  given 


by : 


Qin  = N*(QIO-Kq*PBBAR) 


(H-12) 

Where  N is  the  number  of  fans  supplying  air  to  the 
plenum,  QIO  is  the  equilibrium  air  flow  rate  for  the  fan 
with  zero  gauge  pressure  and  Kg  is  the  slope  of  the 
fans'  airflow  vs  gauge  pressure  curve  . 

Qout  is  airflow  out  of  the  plenum  and  is 
determined  oy  tne  seal  Leakage  eguation: 

Qout  = 2n*Al*  (2*P53AR/?.HOA)  **  1/2 


(H-13) 

Where  Cn  = seal  leakage  coefficient,  Al  = seaL  leakage 
area,  and  RHOA  = mass  density  of  air.  Therefore  mass 
flow  rate  is: 


MBDOT  = RHOA* (Qin-Qout) 


(11-14) 

The  plenum  air  dynamics  are  considered  in 
equilibrium  when  Qin  = 2out 

Aosolute  Plenum  Pressure  is  determined  from  the 
adiabatic  gas  law  to  be: 

?b  = Pa* (MB/  (RH0A*PB3A2)  ) **GAMHA 
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(11-15) 


and , 


P33AH  = ?b-Pa 


( II-  15a) 


Where  Pa  = ambient  atmospheric  pressure  (L3P/F"**2  GAMMA 
= adiaoatic  process  coefficient,  Pb  = plenum  absolute 
pressure,  and  PBBAS  = plenum  gauge  pressure. 


2 . Buoyant  Forces 


Buoyant  force  is  calculated  by  the  submerged 
volume  multiplied  by  the  mass  density  of  the  fluid  times 
the  gravitational  constant. 

Therefore  forward  ouoyant  force  is: 

HBF  = -2*HHO*G* V OL 

(11-16) 

Where  volume  is  given  by  egn.  II-2,  330  is  the 
mass  density  of  watec  (slugs/ft**3)  and  G is  tae 
gravitational  constant.  The  factor  of  2 is  introduced 
to  account  for  2 sidewalls.  The  total  expression  for 
HBF  is: 

HBF  = -2*RHO*3*L1*LD3AR1 * (LDBAR1 / (2*T AM  ( DR  1 ) ) fWSOI) 

(11-17) 

Following  the  same  pattern. 
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ri  3 A = -2*R HO*G*L 2*LDBA 32*  (LD3A  32/ (2*TAN  (D32)  ) =<<510) 

(11-18) 

3.  Beal  Facets 

Forward  seal  force  is  determined  from  egn.  II-9  to 

be : 

Hi  F = -P33Aa*»I3TH*  (LD-L3*TAN  (THETA)  ) /SIN  (31) 

% 

(11-19) 

where  HSF  is  tae  nea  ve  force  due  to  forward  seal. 

In  a siailar  aanner  HSA  is: 

HS A = (?3BAS*?DI?F)  (LD*i,4*TAS  (THETA) /5IN  (32) 

(11-20) 


4 . Planing 

Hclntyre  introduced  the  planing  force  phenomena  in 
his  aodel  to  account  for  the  proven  fact  that  the  CAB 
SES  tends  to  decrease  it's  draft  and  pitch  down  as  speed 
increases.  His  development  of  the  planing  equation  was 
not  well  explained  aid  the  basis  for  selection  of  the 
planing  coefficient  was  not  clearly  stated.  Seidel  Cl  2] 
did  a special  investigation  of  the  planing  force 
phenomena.  His  results  were  based  on  flat  plate  pianino 
effects  to  approximate  the  fiat  bottom  of  the  craft 
s idewails . 
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Since  the  water  is  in  contact  with  the  vertical 
portions  of  the  sidewalls  it  was  felt  a more  accurate 
representation  was  neeced.  Barnaoy  Cl  d]  states  that 
water  contact  on  the  vertical  sides  of  a laning  hull 
rapidly  reduces  planing  lift  advantage  since  one  of  the 
advantages  of  a planing  hull  is  reduced  drag  due  to 
wetted  surface.  planing  Lift  results  in  a speed 
sensitive  tradeoff  between  static  lift  (buoyant  effects) 
and  dynamic  lift  (planing  effect)  that  is  directly 
proportional  to  tne  angle  of  attacn  (for  small  planing 
angles)  and  tne  sguare  of  the  velocity.  Therefore  the 
planing  force  is  of  tne  form: 

HPLAN  = K * V* * 2 *TH  E T A 

i 

From  Seidel's  work, 

I 

HPLAN  = J . 5*AREA*V**2*PI*SIN (THETA) 

(ii-:  i) 

I 

for  saall  angla  approxiaa tions  tais  siaplifias  ro: 

I 

HPLAN  = RH0*A3EA*V**2*PI*TH2TA 

Which  is  the  planing  force  for  a fiat  plate.  3am abv 
gives  no  simple  formula  for  the  losses  due  to  immersion 
so  we  are  left  with  an  arbitrary  loss  coefficient, 

PLCOEF  such  that: 

» 

HPLAN  " (PLCOEF)  *2*HHO*.ARSA*V**2*PI*rHETA 

(11-22) 

The  added  factor  of  '2  accounts  for  2 sidewalls. 

The  planing  area  is  taken  from  the  scale  drawings 


25 


developed  from  sidewall  data,  and  the  force  is  assumed 
to  act  at  the  centroid  of  the  flat  keel  surface,  which 
is  consistant  with  Seidel’s  work.  The  loss  coefficient 
will  be  developed  in  the  section  on  achieving 
equilibrium. 


5.  PllBua  Gauge  Pressure  Lift 


Combining  equations  Il-lla  ,II-11b  , and  II-llc 

results  in: 

HP8ES  = - P 3BAK * W IDTH * (L? W L- (LPWL-LPR) * (LD/BUBHGT)  ) 

(H-23) 

6 . Total  Lift  For ces 

Lift  Total  = HPRSS  * H3F  ♦ H3  A ♦ HS?  +•  HS A + HPLAN 

whicn  must  exactly  equal  craft  weight  (W)  for  the  hull 
to  be  in  vertical  equilibrium,  or; 

tf  <•  HPRSS  + HEF  *■  H3  A «•  HSF  HS  A +■  HPLAN  = 0 

(11-24) 

defines  equilibrium  condition  in  heave  motion. 

H.  CRAFT  PITCH  MOMENT  EQUATIONS 

t . Buoyant  Moments 


1 

f 


: 


14 


I 
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Forward  buoyant  moment  is  the  forward  buoyant  lift 
tiaes  the  effective  lever  am,  or; 


P3F  = - HBF*  (Lever  arm  forward) 


PBF  = -HBF*L5 


(11-25) 


Note  that  the  negative  sign  results  from  the 
definition  of  positive  noments  and  the  fact  that  upward 
forces  are  negative  by  definition. 

P3A  = HBA*LEV2R  AS  ii 


P 3 A = hBA  * L6 


(11-2  6) 


using  the  same  reasoning,  after  buoyant  moment  is; 


? B total  = PBF  * ?3 A 

(II-27) 

2 . Seal  Mg  men  ts 

Since  the  seals  are  relatively  distant  from  the 
C.d.  dynamic  variations  in  the  lever  arms  are  considered 
in  seal  moment  equations  and  are  developed  nere. 


i 


c 


i- 
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1 


I 


The  seaL  force  is  modeled  to  act  at  the 
longitudinal  canter  ofthe  seal  wetted  length.  Therefore 
the  effective  lever  arn  tor  the  seal  moment  is: 

PLSP  * L3  - XSEAL1/2 

where  PLSF  = moment  lever  arm  for  seal  forward.  As  a 
result : 

PS  F = -HSF*PLSF 

PS F = -HSF*  (L3-  (XSEAL1/2)  ) 


(11-28) 


In  a simiLar  manner 


PLSA  = L4  ♦ XSEAL2/2 


PSA  = HSA  * PLS  A 


PSA  = riSA*  (L4*XSEAL2  /2) 


( II  “24) 


PS  TOTAL  = PSF«-?SA 

Note  that  the  dynamic  pitch  sensitivity  te  moment 
arm  length  is  included  in  the  XS SAL  terms  as  lefined  in 
equations  II-S  and  II-3. 
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3 • P 1-ini lowest 

For  the  planing  force  devalopad  in  equation  11-22 
the  planing  force  is  concentrate!  at  the  centroid  of  the 
flat  iceel  which  results  in: 

PPL AN  = L7*HPLAN 

(11-30) 


4 . £i J0.il.  pressur a latent 

First  the  dynamic  variation  in  effective  center  of 
pressure  will  be  developed,  (see  fig. (4)) 

plenufl  length  is  tasen  from  equation  Il-lla  to  be 

LP  = LPML- (LPWL-LP3)  *(LJ/9UBH5T) 

Defining  LDIFF  to  be 

LDIFF  = LPML-LP5 


(ii-3 1) 


and. 


XCP  = XCP0  + (LDIFFOLD/fl'JBHGr)  /2 


(H-32) 


where 


29 


■ ■ ■■  ■«.«  *•+** 


XCPO  = 1PWL/2-XCG 


(H-33) 

and  XCG  is  the  distance  from  the  transom  to  the  C.G. 
(see  fig.  4 .)  Finally, 


PPBES  = -H?RES*XCP 


(11-34) 


5.  Pit £n  Damping.  do  man t 

Any  frictional  damping  terms  dependent  apon  the 
square  of  pitch  rate  (i.e.  tangential  velocity  about  the 
'£  axis)  disappear  when  the  Taylor  Series  expansion  is 
carried  about  the  zero  rate  equilibrium  point.  Since 
for  small  signal  approximation  the  model  should  reflect 
actual  craft  characteristics,  it  was  hoped  that  the 
changes  in  seal  and  plenum  modeling  would  reveal  the 
missing  damping  characteristics.  As  it  turned  out  this 
was  not  the  case.  Examination  of  the  hydrodynamic 
aquations  in  C9J  revealed  an  added  mass  effect  dependant 
upon  a cross  velocity  term  between  linear  X direction 
velocity  U,  pitch  rata  Q,  and  added-mass  at  the  stern 
A335  due  to  slender  body  theory  applied  to  the  sidewalls 
. This  moment  is: 


FP  = -A33S*XLSS**2*U*THErADOT 


(11-35) 

tfhara  XLSS  is  the  distance  from  the  C.G.  to  the  stern, 
and  THETADOT  is  pitch  rate. 
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A3  35  =*  YSS**2*RHO*PI/8 

(11-36) 

where  YSS  is  the  :ross  section  width  of  the  side 
wail  at  the  stern  watecLine.  Therefore, 

FPOANP  = -2*XLSS**2*A33S*U*THETA30T 

= -2*PI*XL3S**2*YLS3**2*8rtO",U*TH2rD3T  /8 

(11-37) 

The  factor  of  2 is  introduced  to  account  for  two 
sidewalls. 

6 . lot^d, 

Since  the  moments  are  all  referenced  to  the  craft 
center  of  gravity,  the  craft  weight  contributes  zero 
moment  in  that  frame  of  reference.  Therefore  total 
moments  are: 

PSP  ♦ PSA  PBF  «■  P3 A PRES  ♦ PPLAN  *■  PDAdP  = 3 

(11-38) 


I.  HEAVE  ACCELERATION  EQUATIONS. 

From  Newton's  law,  acceleration  and  force  are  related 
by: 
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fohcs 


NASS* ACCELERATION 


For  the  linearized  model  we  are  interested  in  the 
acceleration  on  the  craft  and  as  a result: 


ACCELERATION  = FORCS/  MASS  and  therefore  the  heave 
equation  of  notion  becomes: 


ZDDOT  = (HBF*HBA*HSA*HPLAN«-HPRES) /NASS 


(11-39) 

where  ZDDOT  is  the  second  derivative  of  z witi  respect 
to  time. 


J.  PITCH  ACCELERATION  EQUATIONS 


Angular  acceleration  is: 

THETADDOT  = PTOTAL/IYY 

Nhere  ITT  is  the  craft  inertial  moment  about  tie  r axis. 
As  a result  total  angular  acceleration  is; 

THETADDOT  = (PBF  *PB  A +PS  F PS  A*  PPL  AN  *P  PR  ES+PDAMP)  /IYY 

(11-40) 


This  completes  the  derivation  of  all  nonlinear 
aquations  needed  to  describe  tne  t wo- degree-of- freedom 
CAS  5E5.  In  the  following  chapter  these  nonlinear 
equations  will  be  linearized  using  the  Taylor  Series 
expansion  about  a steady  state  operating  point  which 
will  result  in  the  equations  needed  to  generate  a 
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iynaaic  iigitaL-coa^utac  siaalation  of  craft  aation  in 
the  tioe-doaain. 
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III. 


LINSARI2AH21  OF  THE 


A.  TAYLOR  SERIES  EXPANSION 


For  small  perturbations  about  a steady  stats  opera*ing 
point  the  nonlinear  function  2 can  oe  approximated  by: 

Z ~ Z (X,Y) 


Z(0)*dZ  = Z(X(0)  ,Y(0))  * (AZ/iX)  MZ*  (*Z/dY)  *dY 

where  (}Z/bX)  and  (bZ/dX)  are  the  partial  derivatives 
of  tne  function  Z evaluated  at  the  steady  state 
operating  point  X(0),Y(0).  Zancalling  the  steaiy  state 
values  from  Doth  sides  of  the  equation  yields  the 
differential  model  equations: 

dZ  = (bZ/bX)  *dX+-  ( bZ/bX)  * d Y 


and 


Z = Z (0)  *dZ 

where  dZ,dX,  and  dY  ace  the  differential  variables  for 
Z,  X,  Y respectively,  (}Z/6X)  and  (AZ/3Y)  are  the 
partial  derivatives  of  Z evaluated  at  the  steady-state 
operating  point. 
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8.  LIN  EAR IZATION  OF  THE  MODEL  EQUATIONS 


There  are  basically  four  equations  to  linearize  for 
the  t*o-deqree- of- f r eedoi  aodel; 

ODDOT  =*  ? (2, THETA,  M3) 


THETADDOT  = 0 THETA,  rd2TADOT,M3) 


M 300  T = H (2, M3) 


? 33  AH  * 1(2,  MB) 

therefore  the  equations  for  the  differentiiL  state 
variables  are: 


d 20  00  T = (dF/iZ)  «d2«-  (iF/dTHSTA)  *drHSTA<- (iF/iM3)  « JM3 

(III-1) 


iThSTAOOT  * (iO/dZ)  • d2»  (13/dTHETA)  *iTHETA»  (.lO/dM3|  *dM3 

(III-2I 


1M30J  T * ( dH/dZ)  «dZ  *■  (i  1'dTildTA)  MTUST  A*  ('iH/dMBl  MMS 

(XII-3) 
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d?33A3  = (dl/dZ)  *dZ*  (dl/dTHEIA)  *drHETA«-  (dI/dM3|  MM3 

(in-4) 

Note  that  the  left  hand  si  la  variables  are 
differentials,  not  tiae  derivatives. 


C.  LI N EA R I2ZD  EQUATIONS 


1.  Pie :>  ua  Air  Mass  fie*  gee ivatrv as 


Qia  =*  N*(\2iO-K^*P33A  a) 


err- tjj 


Qoat  = Cn*Al*  (2*PB3Aa/3H3A)  **1/2 


(II-1  3) 


N3D0T  = a HO  A*  Uin-tfouti 


( II  - 1 4) 

as  a result  the  differential  ejuatun  for  .13D0T  becomes: 


ddiDOr  = a3CA»  (d^2in- djeat ) 


(III-5) 


l^m  = (djin/dl)  *dZ*  (dQLa/d?B9AR)  • IP 3BAR ♦ndOia/ 11 3 ) *dN3 


(HI-6) 


dQi n/d  P33  A H = -N*Kq 

The  partial  derivatives  ot  Qin  with  respect  to  all  other 
variables  are  zero. 


djout  = 

(d^out/iZ)  *dZ«-djout/dP38AH)  *dPBBA3*  (dJout/ifJB)  M.'IB 


dJout/dP33 A3  = Zr*M*  ( (3  30A/2*  PB3A3)  **  1 / 2)  /RHOA 

(III-7) 

The  partial  derivatives  of  ^out  * itn  respect  to  all 
other  variables  are  zero.  Inserting  equations  EZI-6,and 
III-7  into  II- 14  yields  the  equation  for  the 
differential  variable  d.1 3 DOT : 


ih3D0r  = (-8H[3A*N*Kq-Ca*U*  ( (SHOA/2*  ?S3  A 3)  *«1/2)  ) «dPB9A3 

(11-15) 


n ua  gauge  pressure  Derivatives^ 


Pb  = ?a*(ab/3HDA*Vb) 


and 


(rn-S) 


P d a A 3 


?b-?a 


Since  Pa  is  considered  constant  witnin  the  time  frame  of 
this  model. 


dP33AR  * dPb 


(IIT-9) 


As  a resalt, 

d?  3 3A  3 = 

SAHA*  Pa*  (.13  (0)  /3H0A*Vb  (3)  ) **  (3  All  A-  1)  *d  (1B/RH0A*VB) 
where 

d (13/ShOA  *V3)  = ( 1 /R HOA <■  V 3 ( ))  ) * d 13+ ( 13/3  dO  A)  « i ( 1 /V  B) 

d (13/RHOA*  Vi>)  = (1/3HOA*V3(0)  ) * d (1 3-  ( 1 3 / 3 HO  A * V 3 * * 2 ) * d V 3 


dV3  = d (V  N- A3*  LD) 

dV3  = -A3  (0)  “dlD/dZ-LD  (d)  *dA3/iZ 
remember! aq  that; 

LD  = Z + Z5 
then 

dLD/dZ  = 1 


we  now  need  the  partial  of  effective 
area  wit a respect  to  2.  Defining 


plenum  roof 


A3 


AoWi- ( A3WL- A3R)  ‘L0/5UBH3T 


( 1 1 1 - 1 0) 

A3  41.  = WIDTH*  LPW  L 

A33  = WIDTH*L?R 
further  definin j 


A3JIFF  = A3WL-A33 

which  is  a constant,  tie  effective  plenum  roof  area  is: 


A3  = A5WL- A5DIF?«LD/5U3H3r 

usin*  the  previous  definitions  trie  partial  of  affective 
roof  area  becomes: 


d A3/i u = -A3DIFF/3U3H3P 


and  tne  differential  voLusne  is: 


d Vb  =-  (A3  (0)  -LD  (0)  *2  * (A30 1 FF/3U3HviT)  ) • iZ 


with  the  result  that 


d (,13/riHOA*  Vb)  = (1/HHOA*Vb  ( 0)  ) *dd  «• 


(HI-  1 1) 


( I II-  1 2) 


(ho  (0)  /3HOA*Vb  (0)  **2)  * (AB  (0)  -2*LD  (0)  * (A3DIFF/BJ3HGT)  ) *dZ 


3 b 


t 

f 

h 

H 

n 


f 


x 


1 


fia ally , 


II 


dPBBAH  = GAHMA*Pb(3)  * ( (1  / Mb  (0)  ) *d!10  *■ 

(1/'/b(0)  * (AB  (0)  -2*LD  (0)  * (A  B DI FF  /3  U3HGT)  ) ) *dZ 

(III-1  3) 


3.  Heave  Force  Derivatives 


The  total  heave  force  equation  is: 

H (total)  = H5F+HBA+HSF+H5A+HP1AN+HPRES 

The  total  partial  derivatives  will  be  done  term  by 

term. 

For  the  forward  buoyant  force  the  following 
simplifications  are  made: 


H3F  = -2*RH0*G*L1*LDBAR1*  (LDB AR1 / (2* TAN  ( DH  1 ) ) a- W 3 1 0 ) 
let  K 1 = 2*RHO*G*L1  then, 

HBP  = -K1  * LD3AR  7*(LDBAH1/(2*TAN  (DR  1 1 ) + US  1 0) 

where  LDBAR1  = (LD- L5TA N (THETA) ) . Combining  the 

previous  two  equations, 

HB?  = -K 1 * ( LD-15  *TA  N (THETA)  ) * 

( (L  D-L5  *TAM  (THETA)  / (2  * TAM  (DR1)  ) «•  rfS  1 0) 


assuming  small  pitra  angles,  TAN  (THETA)  can  be 
approximated  by  THETA  which  results  in  the  following 
simplir ication : 


HBF  = 

-K1*(LD-L5*THErA)*((La“L5*THETA)/(2*rAU(DR1))«-W310) 

Talcing  tne  partials  of  HBP  with  respect  to  Z and 
THETA  yields: 

i (HBF)  /id  = -K1*  ({ID  (0) -L5*TH2TA  i0)  /TAM  (Dai)  ♦WS  10 

(ITT-14) 


and 


d (HBF)  /dTHETA  = 

E 1*L5*  ( (LD  (0)  -L 5* THETA  (0)  ) / (2* TA  N (D/r  1 ) ) «-WS10) 
♦ K1*  (-L5/(2*TAN  (DS1M  * (LD(0)  -L5*  THETA  (0)  ) 

Which  simplifies  to  : 


d (HBF) /dTHETA  = 

K 1 * L5  * ( ( L D ( 0 ) -L5  *Tf!ETA  ( 0)  ) / (T  AM  (DR  1 ) ) ♦ 4 3 1 0) 

(I  II-  15) 

Using  the  same  aporoach,  the  partials  for  buoyancy 
force  aft  are: 

d(H3A)/dZ  = -K2*  ((LD  (0)  *-L6*THETA  (l))  ) /(TAN  (DK2))  *WS20) 


4 1 


( 1 1 1 - 1 1' ) 


d (HBF)  /dXHETA  = 

-Lo«>K2*((LD(O)*Lb*rHETA(O)/(TAfl(0R2))«-«IS1O) 


(ITT  — 17) 

For  the  forward  saal  the  total  differentials  are; 

DHSF  =(d(HSF)  /dZ)  *dZ ♦ (d  (H  SF)  /d  TH  E I A)  *d  TH  ET 
«■  (d  (HSF)  /dPBBAfi)  MP3BAR 

the  following  simplification  is  made; 

dSF  = -P3flAa*aiDTH* (LD-L3*TAM (THETA)  /SIM (3 1) 

(11-19) 

let  HIDTH/SIM (31)  = K3  and  TAM (THETA)  = THETA,  then, 

H5F  = -K3*?3BAR* (LD-L3* THETA) 
with  the  resultant  partials  being: 

d(KSF)/dZ  = -K3*PBBAR  (0) 

(III-18) 


i (HSF) /dTHSTA  = K3*L3*P8B  AR  (0) 


(IIT-19) 


i (HSF)  /dPBBAR  = -K3*  (LD  (0)  -L3*THETA  (0)  ) 


(III-20) 
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The  partials  of  after  seal  farce  are  founl  in  the 
same  manner  to  be: 

dHSA  = (d  (HSA)  /dZ)  *dZ*  (d  (HSA)  /dTHSTA)  *dTHETA 
*•  ( d (HS  A/dP  BB  AR)  * dP  B 3 A S 


HSA  = - (?BBAR*PDIFF)  *WrOTH*  (LD*L4«TAN  (THETA)/SIN  (32) 

(11-20) 

let  K 4 = W IDTH/SIN (3  2)  and  assume  small  TH ST A so  that; 


HSA  = -K4* (PBBAStPDIPF) « (LP+L4*TH3TA) 
the  partials  then  become: 

d (HSA)  /dZ  = -Kh*  (PB3AR  (3)  *PDTFF) 


(III-21) 


d (HSA) /dTHSTA  = -K4*  (PB  A 3 (0)  «■  PDIFP)  * L4 


(1 11-22) 


a (HSA)  /dP B BA 3 = - K4 * (LD ( 0)  *L4 *TH ETA (0) ) 


(1 11-23) 


Next  the  planing  force  partials  are  derived: 


H PL  AN  = -2*PL20EF*RH0*'USA*V**2*PI*rH2rA 
let  KPLAM  = (2*?LCOE?F*RHO*AREA*PI)  then. 
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1 


HPLAN  = -KPLAN  *V**2*THETA 
and , 


d (HPLAN)  = (d  (HPLAN)  /dTHET  A)  * dTHET A 
for  this  constant  speed  model,  therefore, 

d (HPLAN) /dTHETA  = -KPLAN  *V  (0)  **2 

Plenum  Gauge  Pressure  lift  w as  defined  to  be; 


( r ii-24) 


HPRE5  = - AB*PB6AR 


as  a result; 


d (HPRES)  /dPBSAR  =■  -A3(3) 


(II-  11c) 


(XIX— 25) 


4 • Pitch  Moment  Partials 


The  total  moments  acting  on  the  craft  are: 


PTOTAL  = PBF*PBA*PSF+P3A*PPLAN+?PR2S*PDAMP 


(H-38) 


The  partials  of  11-38  will  be  found  term  oy  term, 
Deginaing  with  the  forward  buoyant  moment: 
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.. . 


PBF 


-HoF*L5 


therefore: 


d (?BF)  =»  -L5*d  ( HB F) 
so  that: 


(11-23) 


d(P3F)  = -L5*  (d  (HBF)  /iZ)  *dZ-L5*  ( (i  (H3F)  /dTHETA)  MTHETA 

(III-2b) 

and  the  partial  of  the  after  buoyant  moment  is: 

P3A  = Lt>*HBA 


(II-2b) 


dPSA  = Lo*  dHBA 
tneraf  ore , 

d(?BA)  = L6*  (d  (HBA)  / iZ)  * dZ»L6*  (d  (H3A)  /dTHETA)  *d THETA 

(111-27) 


Seal'  forces  partials  must  include  the  dynamic 
variation  in  lever  arm  length,  waich  for  the  forward 
seal  is: 

PSP  = - HSF*PLS  F 
or, 

I 


PSF 


(P3BAR  *A5EAL1 ) *PL5F 


therefore , 

d (PSF ) /dP3  BAR  = ASEAL1  (0)  *PLSF  (0) 

(III-29) 

and 

d (PSP)  dZ  = PBBA3  (0)  *d  (?LSF*ASEAL1)  /iZ 
where 

d (ASEAL  1*PLSF)  /dZ  = PLSP  (0)  *d  (ASEALl ) /dZ  ♦ 

ASEAL 1 (0)  *d (PLSP) /iZ 

AS  SAL  1 = WIDTH* ( LD-L3*T 3 ET  A) /SIN (3 1) 
and 

PL3F  = L3- (LD-L3*THETA)  /2*SIN  (31) 
therefore , 

d(?LS?)/iZ  = - 1/  (2*3  IN  (3  1) 
and 

I(A52AL1)/dZ  = WIDTH /SIN  (31) 
as  a result, 

d (PSP)  /dZ  = PB3AR  (0)  *(  (PLSF  (0)  *WIDTH) /SIN  (3  1)  - 
ASEAL  1 (0)  / ( 2*SIN  (31)  ) 


U6 


( 1 11-29) 


d (psf)  /dTHETA  = PBAR  (0|  *d(PLSF*ASEAL1)  /dTHETA 

d (PLSF*ASEAL1)  /dTHETA  = 

PLSF  (0)  *d  (AS5AL  1 ) /dTHETA  ASZAL  1 (0)  *d  (PL S?)  /dTHETA 

d (AS2AL1)  /dTHETA  = - L3*rf I DTH/S IN  ( 3 1 ) 

d (PLSF) /dTHETA  = L3  / (2  * S I N (31) 
as  a result, 

d (PSF)  /dTHETA  = -P33AR  (0 ) * ( (PLSF  (3)  *L3*VID?H)  SSLS  (3  1)- 
ASEAL1  (0)  *L3/  (2*SIN  (3  1 ) ) 

(I II -30) 

For  the  after  seal  the  derivatives  proceei  in  the 
manner: 

PSA  = - (PB3AH  ♦ PDIFF)  *A5SAL2*PLSA 


d (PSA) /dPEBAR  = -ASE AL2 (0) *PISA  (0) 


(III-31) 

In  the  same  manner  as  shown  in  the  previous 
paragraph,  the  partials  of  moment  due  to  seals  aft  with 
respect  to  the  other  variables  are: 


h7 


— 


d (psa)  /dZ  = -P38 AR  (0)  * ( ( PLSA  (0)  /SIN  (32)  * 

AS  EAL2  (0)  / ( 2*SI N (32)  ) 


( I II -32) 


d (PSA) /dTHETA  = -PBAR(3l  * ( (PLS  A(0)  *L4*MIDTH) /SIM)  32)  ♦ 
ASEAL2  (0)  *L4/  ( 2*S  IN  (3  2)  ) 

( 1 1 1 -3  3) 


The  partial  of  planing  moment  is: 


PPLAN  = L7*HPLAN 


(11-30) 


d(PPLAN)  » L7*d  v riPLAN) 


i (PPLAN) /dTHETA  = -L7»KPL AS*V ( J)  *»2 


( 1 1 1 - 3 4 ) 


Par  this  constant  speed  model  all  otner  oartials  are 
identically  zero.  The  Plenum  dauge  pressure  also  has 
dynamic  variation  in  the  moment  lever  ara  and  the 
derivatives  are: 


PPRSS  = - XCP*  H P R SS 


the  ref  ore , 


(11-34) 


dPPRSS  = -HPR2S  (0) d (XCP)  -XCP  (0) *i (HPRE3) 


¥ nere 


XC?  = XC?0>  (LDIFP*L3)  / (2«3U3Hijr) 
t heref oci, 


d (XCP)  /dZ  = LDI?F/(2*B'J3HGT) 
a n J 


d(?PR2S)  = - ( (HPSES  (0)  *LDIFP)  / (2*BUBHGr)  ) *d2 
-XCPlO) *d (HPRES) 

substituting  equation  (EII-25)  yields: 


d (?  PR  2 5)  = -(  (H?RE3(0>  «LDIFF)  /(2«3JBHGD  ) 
♦XCP(O) * A3 (3) *dP3oAR 


Finally,  the  partial  of  pitch  damping  monant  is: 


PDAM?  = -2*PI*XLSS**2*YL5S**2*RH3*'.l*  THETA  DOT/9 


(TI-37) 


let 


2*Pr*XLSS**2*YSS**2»RHO/3  = KDAflP 
then 

PDAhP  = -KDAMP* V*THETADDT 

as  a result,  since  this  is  a constant  speed  aodal. 
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d l ? 0 A h F ) / iTHSTOr  = -KDA'1?*V  (0) 

(HI-35) 

wh-»re  THETDr  is  the  first  derivative  of  theta  with 
respect  to  time. 

5.  £ot&l  2i t f er jati 1 1 3 

Since  ? a h • A and  tae  interest  is  in 
acceleration  terns,  A = F/.1  and  dA  = i(?)/S.  tne  result 
beuj  that  all  differentials  must  oe  divided  by  craft 
mass  to  acaieve  the  desired  linearized  heave  results. 
The  same  general  comment  holds  true  for  moment 
differentials,  that  is;  they  all  must  be  divided  by  the 
craft  moment  of  inertia  about  the  Y-axis. 

All  differentials  must  be  sorted  by  the  variable 
of  differentiation  so  that  the  linearized  differential 
equations  of  motion  may  be  written.  As  an  example; 

iZDDJT  = { (d  (HB  F)  /id)  ♦ ( i { H 3 A ) /id)  ) /MASS)  *1Z  + 


The  State  Variable  representation 
will  be  presented  in  a later  section 
however  it  is  convenient  to  define  the 
at  this  tiae. 


of  the  system 
of  this  chapter. 
State  Variables 


t 

X = (iS#dSCOr,drHETA  , irrlSTDT,  ih3) 


( 111-30 

and 


SO 


I * 1PB3AR 


( r II- 37) 

since  P9BAR  is  tint  a state  viruoLe  tne  definition  for 
dPSBAR  aast  be  substituted  for  d?33AR  in  all  equations. 

o.  Tot  al  deave  Force  different  ills  in  2 


Here  the  total  differentials  will  be  found.  No 
explanatory  text  will  be  included  unless  special 
saDstitut ions  requiring  an  aside  are  necessary. 

The  cuoyancy  terms  are: 

3H32  = d ( H3F)  /dZ  ♦ d(H9\)/dZ 


(111-38) 

A term  oust  be  included  which  was  not  previously 
derived  to  account  for  the  sensitivity  of  seal  forces 
due  to  P9BAR  since  it  is  not  a state  Variable. 

DHSPZ  = (d  (HSF)  /dP9BAR)  * (d  (PB8AR)  /iZ) 
as  a result, 

d (HSF) /dP33AR  = -K 3*  (LO  ( 0 ) -L3 *TK  ET  A ( 0)  ) 


(111-20) 


and 


DHSFPZ  = -K3«  (LD  (0)  -L3*THSTA(  0)  ) «d  (P3BAR) /d2 


siailarly  for  the  after  seal  moment: 


i (HSA)  /dP  3BAH  =-  K4*  ( LO  (3  I ♦L4*THETA  (0)  ) 


(1 11-23) 


and 

OSSAPZ  = - K4* (LD  (0)  +L4*P33PA(0))  *i(?B3AR)/dZ 
rha  total  seal  force  differential  in  Z becomes: 

DtiSZ  = l(HSF)/dZ  + d(H3A)/dZ  OHSFPZ  * DHSAPZ 

(III-39) 


For  plenum  gauge  pressure  lift  differential  the 
sane  substitution  must  be  made. 

dd?32S  = - A 3 ( 0) * dP  33  AR 
therefore , 

3HPZ  = - A B (3)  * d ( P S B A 3)  / i Z 

(rrr-40) 

Define  the  sensitivity  coefficient  of  ZDD'JT  to  Z 

to  oe : 

DZZ  = (DHBZ  ♦ D3SZ  ♦ DHPZ) /SA5S 

(III  -4  1 ) 


Total  Dif f eraatials  in  THETA 


7. 


The  following  partials  are  now  defined: 

DH3IH  = d (HBF) /ITHETA  ♦ d (HBA)  /dTHSTA 

(III-42) 

DH3TH  = d (HS?)  /dTHETA  +i  (H3A)  /dTKETA 

(IIT-43) 

DHPTH  = d (HPLAN)  /dTHETA 

(III-44) 

Define  the  sensitivity  coefficient  of  2DDDT  due  to 
THZTA  to  be: 

DSTH  = (DHBTH  * DHSTH  ♦ DHPTH)  /WAiS 

(III-45) 


3.  Total  Differentials  due  to  13 

Here  since  dP33A3  is  not  a state  variable  the 
partial  of  PBBAR  with  respect  to  .18  must  be  suostituted 
in  the  plenum  lift  and  seals  lift  terms. 


d (HS?) /dPBBAR  = *K3*  (LD  (0)  -L3* THZTA  (0)  ) 


(III  — 20) 


oy  substituting  d(PB3A3l  /d13  for  1P33AR  the  result  is: 


f 


53 

1 ' Ml 


I 


d ( H S F P ) /dMB  = -K3*  (LD  (J»  -L3*THETA (0)  ) *i (PB3AR)  /id3 


(iir-4b) 


similarly , 


d (HSAP)  /dMB  = -K4*  (LD(0»  f L4*THETA(0|  ) *d  (PB3AR)  /dMB 


( 1 1 1 -4  7) 


d ( dPSES)  = -A3  (0)  *dp 3BAR 


usin^  the  previous  substitution  yields: 


d (HP8E5)  /dMB  = -A3  (0 ) * d (P  3 BAS)  /dMB 


(1 11-48) 


The  sensitivity  coefficient  of  heave  due  to  MB  can 
now  be  defined  to  be: 


Dd.13  = (d  (HSFP)  / dMB  +d  (HSAP)  / i .13  * d (HPSES) /dMB)  /MASS 


(I  11-4  4) 


3 • X2iLll  Pitch  Moment  Differentials  in  Z 


DP3Z  = d (F  BF)  /dZ  +■  d (PSA)  /dZ 

(I  II -50) 

The  dPBBAR  substitution  yields  the  partial  of 
seals  with  respect  to  pressure  wita  respect  to  Z. 


DPSFPZ  = ASEAL1  (0)  *PLSF  (0)  *d  (PBBAIt)  /dZ 


(I II -5 1 > 


DPSAPZ  = -ASEAL2  (0)  *PLSA  (0)  *d  (PBBAR)  /dZ 


( 1 1 1 -5  2) 


therefore  the  total  seal  differential  in  Z becones: 


DP SZ  = DPSFPZ  ♦ DPSAPZ  fd(P3F)/dZ  + l(?SA)/dZ 


(III  *5  3) 


Next  the  Plenum  Zuge  pressure  term  is  derived, 


d (PPaES) /dZ  = -HPRES  (U|  * L D IFF / ( 2 * BU 5 HGT ) 


d (PPRES)  /d?3DA3  = AB  (0)  * YCP  (0) 


theretore. 


DPPRESZ  = A3  (0)  *XCP(0)  «d  (PBBAR)  /dZ 


DPR  ES  Z = d(?PRSS)/dZ  + d(  P PRES  Z) /JZ 


(111-54) 


The  sensitivity  coefficient  of  THETA  wit  a respect 


to  Z becomes: 


DT  HZ  = (DP5Z  *■  DPSZ  * DPPRES)/IYY 


( 1 1 1 -5  5) 


10.  Total  floa  en  t Diffecentiels  in  THETA 

DP3TH  = (d  (P3F) /dTHETA  ♦ d (PB  A ) /dTH  3 TA) 

(III-56) 

DP5TH  = (d  (PSF)  /dTHETA  «•  1 (PSA)  /iTHE  TA 

( r ii  ~57) 

DPPLTH  = d (PP  LA  N ) /dTHETA 

( 1 1 1-5  8 ) 

define  the  sensitivity  coefficient  for  THETA  to  be: 

DTHTH  = (DPSTH  *■  DPSTH  *■  DPPLTH)  /I'll 

(1 11-5  9) 

Also,  define  the  sensitivity  coefficient  in  THETDT 

to  be : 

DTHDTH  = (d  (PDAHP)  /dTHDT)  /IYY 

(I  II -60) 

11.  To  t.i  1 doaent  Pit  fere  nr  ills  in  d3 


5o 


d P3 AM3  = (d  (PSA)  /dPBBAR)  * (d  (PBBA3)  /dMB) 


( 1 1 1 -6  1 ) 


dP3FMB  =(d  (PSP)  /1PB3AR)  * (d  (P33AR)  /id 3) 


( 1 1 1 -6  2) 


QP5M3  = DPSAMB  * DPSFMB 


(III— 63) 

Inserting  the  dPB3A3  suostitution  yields  the  partial  of 
plenum  gauge  pressure  moment  with  respect  to  M3: 

i ( ?P3 £3) /d P B3AR  = - XCP  (3 ) *DH?R SS 

d (HP3SS) /ddB  = -A3  (0)  *i  (PB3AR)  /d.13 


DPPMB  = XCP(0)*A8(0)  *d  (? 3B  AR)  /dMB 


(III-64) 


Define  the  sensitivity  coefficient  of  THETA  due  to 
M3  to  be: 


DTHMB  = (DPSMB  *■  DPP  MB)  /I  YY 


(TTT-65) 


To t a 1 Differentials  of  Air  Mass 


12. 


d (MB) /dMB  = (d  (MB)  /dPBBAR)  * (d  (PBBAR)  /dMB) 
define  the  sensitivity  coefficient  to  be: 


DM  3d  3 = d (MB) /dMB 


(III-66) 


define  the  sensitivity  coefficient  in  Z to  be: 


DM  BZ  = i(M3)/dZ 


(III-67) 


D.  TOTAL  MODEL  SENSITIVITY  COEFFICIENTS 


1 * Heave  S en  s it  jv  it  y Coef  f icients 


The  sensitivity  coefficients  for  the  Linearized 
differential  equation  in  heave  acceleration  are: 

DZZ , DZTH,  and  DZ  MB. 


2.  Pitch  Sensitivity  Coef f icients 

The  sensitivity  coefficients  for  the  linearized 
differential  equation  in  pitch  acceleration  are: 

DTHZ,  DTHIK,  DTHDTH,  and  DTHM3. 

2*  P le  nam  Air  Mass  Sensitivity  Coefficients 
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The  sensitivity  coefficients  for  plenum  air  mass 

are: 


DMBZ,  and  DM3MB. 

4.  Pjen urn  Gau  qe  Pressure  Sensitivity  Coefficients 

. The  sensitivity  coefficients  for  plenum  Gauge 
Pressure  are: 

DP  3 Z , and  DPS  SB. 


S.  STATS  SPACE  REPRESENTATION 


The  state  matrix  form  of  a dynamic  system  is: 


XDOT  = AX  BU 

where  X = a vector  of  variables  representing  the  states 
of  a set  of  first  order  linear  differential  equations, 
and  J is  the  vector  of  forcing  functions,  and  A and  B 
are  the  matrices  of  the  respective  gains.  The  form  of 
the  equations  is: 


XI  DOT  = X2 


X2D0T  = X 3 


X3DOT  = K 1 *X  1 K2*X  2 «■  K3*X3 


59 


tor  the  linearized  aodel  of  the  CAB  S2S  the  state 
variables  were  defined  to  be  : 

t 

X = (dZ,  dZDOT,  dTHETA,  dTHETDT,  1.10) 


As  a result  the  state  uatrix  representation  of  the 


linearized 

unf or  ced 

Hea/  e 

-Pit  ch 

aodel  is: 

x inor 

0 

1 

0 

0 

0 

XI 

X2U0T 

DZZ 

0 

DZT  H 

0 

DZM3 

X 2 

X3D0T 

= 

0 

0 

0 

1 

0 

XJ 

X4  JOI 

DIHZ 

0 

or  HTH 

3THDTH 

DTHM3 

X4 

x5oor 

DN3Z 

_ 

0 

0 

0 

DM313 

X5 

(III  -6  3) 

and 


dPBBAi  = 0PBZ*X1  * D PB  M B * X 5 
1II-o9 


t»  0 


IV. 


JJTIE  2QMAIN  VALIDATION 


The  linearized  model  developed  in  chapter  III  was 
validated  against  the  the  Six-degree-of- freedom 
nonlinear  model.  There  were  several  changes  made  in  the 
simulation  techniques  used,  and  at  the  same  time  the 
changes  in  center  of  pressure  dynamic  variation  and 
planing  forces  were  entered  into  tne  6 DDF  model. 
Establishment  of  equilibrium  for  the  linear  model  was 
accomplished  off-line  using  a TI-59  programmable 
hand-held  calculator  with  printer. 

The  first  change  in  simulation  technique  is  the 
operating  point  at  whica  the  linear  sensitivity 
coefficients  were  evaluated.  In  dclntyre's  model  these 
coefficients  were  evaluated  at  the  initial  operating 
point.  In  this  work  the  coefficients  are  evaluated  at 
the  final  steady-state  equilibrium  point  since  this  is 
customary  in  linearization  wor'c. 

The  second  major  change  in  simulation  was  the  method 
used  to  determine  the  planing  force  loss  coefficient 
(PLCOEFF)  and  the  establishment  of  the  arbitrary 
location  of  tne  initial  plenum  center  of  pressure.  In 
this  model  all  heave  forces  are  well  defined  (as  modelei 
here)  with  the  exception  of  pianing  force  lift.  As  a 
result,  craft  weight  and  ail  heave  forces  except  that 
due  to  planing  lift  are  summed  and  the  residual  assumed 
to  oe  due  to  the  planing  phemonema.  This  residual  is 
then  used  to  determine  (off-line)  the  value  of  PLCOZFF. 
Note  that  this  method  also  accounts  for  any  lift  due  to 
aerodynamic  forces  on  the  bow  that  were  disregarded. 


M <3 :c t all  moments  are  summed  about  the  C.G.  except  the 
moment  due  to  plenum  pressure  lilt.  The  residual  moment 
is  assumed  to  be  due  to  plenum  pressure.  This  residual 
is  used  (again  off-linet  to  estaolish  the  initial  center 
of  pressure.  This  accounts  roc  any  moments  not 
previously  accounted  for  (such  as  aerodynamic  forces  on 
the  bow.)  McIntyre  used  the  planing  lift  coefficient  to 
force  the  model  into  ajreemeut  with  the  6 DDF  model  for 
natural  pitch  frequency  and  used  a large,  arbitrary,  and 
static  moment  to  force  the  models  into  steady  state 
pi  C w h angle  agreement.  It  was  felt  that  using  the 
present  method  of  leteraination  would  leav^  the 
equations  in  s form  that  would  be  more  amenable  to 
extraction  of  the  factors  most  Iiractly  affecting  craft 
cnaracteristics  as  modeled  herein. 

A.  LINEAR  PROGRAM  DESCRIPTION 


Jser's  instructions  have  neen  included  as  appeniix  \. 
The  program  flow  can  be  isolated  into  four  main 
sections. 

1.  First  the  values  for  craft  weignt,  draft,  pitch 
angle,  plenum  pressure  and  speed  are  entered.  To 
initialize  plenum  air  mass,  draft  is  used  to  sj1v»  for 
initial  plenum  air  volume,  which  is  then  used  with 
initial  plenum  gauge  pressure  to  solve  equation  11-15 
for  plenum  air  mass.  Equation  11-11  is  then  solved  for 
jin.  siout  is  set  equal  to  Qin  and  equation  11-13  is 
solved  for  seal  leakage  area.  All  othei  heave  and  pitch 
nonlinear  equations  are  solved  and  their  values  are 
listed  along  with  cesiluals  of  heave  aid  pitch 
summations.  These  resiluals  are  aided  to  the  initial 
values  of  C.G.  and  angular  acceleration  and  tn e initial 


values  are  saved  for  the  time  domain  simulation  section. 


i 


2.  The  same  operations  are  earned  out  for  tie  steady 
state  values  of  weight,  draft,  piten  angle,  and  plenum 
pressure.  The  residuals,  forces  and  moments  are  again 
listed. 

3.  These  steady  state  values  are  used  to  evaluate 
the  sensitivity  coefficients  developed  in  cnapter  III 
which  are  also  listed. 

4.  Finally  the  state  variable  solution  is  simulated 
using  a four-point  fixed  step  Runge-Kutta  algorithm  and 
listed  for  the  print  interval  desired.  The  option  is 
included  to  output  either  printer  plot  grapnical  display 
using  the  NPS  printer  plot  suoroutine  PLOT?  or  the  X-Y 
Versatec  plot  suoroutine  DRAW?.  The  latter  should  be 
used  for  final  runs  only,  witn  PLOT?  used  on 
intermediate  runs.  Provisions  are  made  for  up  to  300 
discrete  output  points  for  the  time  domain  tabular 
listing . 


3.  5IX-QSGREE-0F-FREEDDM  MODEL  SOLUTION 

Two  changes  were  made  in  the  6 DOF  model  to  account  in 
a more  analytic  manner  for  the  observed  tendency  of  the 
XR-3  to  pitch  down  as  speed  increased. 

The  first  change  was  to  remove  an  empirical  eguatior. 
used  to  correct  the  center  of  pressure  L5]  and  replace 
it  with  the  dynamic  center  of  pressure  variation 
developed  in  chapter  El. 

The  second  change  was  to  include  in  the  sidewall 
subroutine  the  nonlinear  equation  for  planing  forces, 
witn  the  loss  coefficient  adjusted  to  ma<e  up  (as 


! 
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closely  as  possible)  foe  lift  loss  iue  to  the  revised 
modeling  of  the  plenum  chamber. 

Finally,  the  6 DOF  model  was  ran  it  both  o722  lbs  and 
t>050  lbs  weight  to  establish  the  initial  and  final 
operating  points,  followed  by  the  weigat  removal 
transient  run  for  time  domain  validation  of  tie  linear 
heave- pitch  model. 

It  should  be  noted  that  furtasr  study  of  the  6 DO? 
model  is  needed  to  evaluate  the  results  of  the  two 
changes  and  to  more  closely  establish  the  arbitrary 
initial  location  of  the  center  of  pressure.  The  large 
amount  of  CPU  and  turn-around  program  times  re]uired  to 
run  tue  nonlinear  model  precluded  including  tait  study 
in  this  work.  Also,  for  tie  nonlinear  model, 
X-lirection  thrust  was  aeld  constant  during  the  run  with 
the  result  that  speed  was  3d. S3  knots  and  still 
increasing  at  the  end  of  the  five  second  run.  This 
effect  will  be  further  iiscussed  in  the  next  section. 


C.  COMPARISON  OF  SSSOLTS 


The  graphical  output  for  the  five  second  runs  are 
shown  in  figures  5 through  3 ror  tae  linear  model  and 
fioures  9 through  12  for  the  nonlinear  model.  The  fixed 
integration  step  size  was  0.005  sen  for  the  linear  model 
and  0.001  sec  for  the  nonlinear  model. 

1 . Plan ua  dauge  Pressure  Transient 


6U 


(see  fig.  5 and  9> 


In  general  bota  models  oithiaitid  a rapid  initial 
redaction  in  pressure  fallowed  ay  a smooth  transient 
back  towards  the  final  operating  pressure.  For  the 
linear  model  the  minimum  occared  at  100  ailli- seconds, 
and  110  Billi-seconds  for  the  j OOF  model,  with  the 
respective  pressures  being  21.38  and  21.30.  Tie  final 
values  were  24.31  and  24.08  for  the  linear  and  nonlinear 
model  respectively. 

As  can  be  seen  from  the  graphical  data  excellent 
agreement  was  acheived  between  the  two  models. 

2.  Draft  Transient 

(see  fin.  6 and  10) 

Both  models  exhibit  a smooth  exponential  type 
transient  towards  the  final,  lower  steady-state  draft. 
The  linear  model  transient  settles  towards  a lower  value 
(6.04  inches  as  compared  to  6.44)  wnich  is  in  part 
attributed  to  a difference  in  modeling  of  the  seal 
forces.  The  nonlinear  aodel  uses  a complex  combination 
of  hydrodynamic  and  nydrostatic  lift  where  the  linear 
model  uses  a simple  pressure  differential  effect.  Also, 
since  the  nonlinear  model  is  settling  to  a signif icar.tly 
lower  pitch  angle  (see  para.  4 below),  the  planing  force 
effect  is  reduced  in  that  model. 

3.  C .0  Acceleration 

(see  fig.  7 and  11) 

Both  models  show  a rapid  positive  transient  in 
acceleration  that  closely  follows  tne  negative  transien* 
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in  plenum  pressure.  Both  models  reached  the  initial 
peak  positive  C.G.  acceleration  at  approximately  20 
milli-seconds  after  the  plenum  pressure  completed  the 
negative  transient.  The  linear  model  appears  to  have 
greater  coupling  between  pitch  and  heave  looking  at  the 
C.G.  acceleration,  but  this  is  not  actually  supportable 
since  the  linear  model  has  greater  pitch  amplitude 
excursions  about  the  exponential  pitch  transient  to  the 
lower  steady  state  value  (see  fig.  3 and  12.)  In 
general  the  two  models  show  excellent  agreement  in 
transient  and  damped  3.3.  acceleration  behavior. 

4 • Pinch  Transient 

(see  fig.  8 and  12) 

It  is  in  this  area  that  the  two  models  are  in 
greatest  disagreement.  3oth  exhibit  a general  damped 
sinusoidal  behavior  about  a center  exponential  transient 
to  a lower  steady  state  value.  In  both  cases  the 
initial  excursion  is  downward.  The  natural  pitch 
fequency  of  linear  model  is  approximately  4.3  radians 
per  second,  whereas  the  nonlinear  model  exhibits  a 
natural  pitch  frequency  of  6.93  radians  per  second. 
Thera  is  felt  to  be  one  major  factor  contributing  to 
this  large  discrepancy,  which  is  the  inclusion  of  the 
planing  force  in  the  nonlinear  model.  Data  presented  by 
Thaler  and  GerbaD4}  showed  a natural  pitch  freguency  of 
5.3  radians  per  second  and  much  higher  oscillations 
about  the  exponential  transient.  It  is  felt  that  the 
simplified  planing  forces  of  the  linear  model  are 
included  in  the  much  moce  complex  hydrodynamic  equations 
of  the  nonlinear  model  and  are  tnerefore  redundant  in 
that  model.  Due  to  time  constraints,  that  avenue  of 
investigation  is  left  for  future  studies.  The  linear 
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model  is  tending  towards  a steady  state  value  on  the 
order  of  0.45  degrees  whereas  the  nonlinear  model  is 
tending  towards  approximately  0.30  degrees.  This  large 
discrepancy  is  partiaLly  due  to  the  redundant  planing 
force  in  the  nonlinear  model  and  the  fact  that  speed  is 
also  increasing  in  the  nonlinear  model  whicn  tends  to 
accentuate  the  planing  force  and  force  the  craft  to 
pitch  down. 

In  the  next  chapter  a signal  flow  graph  will  be 
developed  and  Mason's  3ain  Rule  applied  to  explain  the 
linear  model  characteristic  cesponse  in  a more 
analytical  fashion. 


o? 


SIGNAL  FLOW  GRAPH  AND  NASON’S  SAIN  HOLE 


V . 


A.  INTRODUCTION 


The  State  Variables  defined  in  Chapt  II  are  used  to 
develop  a signal  flow  graph  of  the  craft  heave-pitch 
dynamics.  Nason's  Gaia  rule  is  applied  to  the  signal 
flow  grapn  to  derive  analytical  5-domain  polynomials 
expressing  the  transfer  functions  for  pitch  moment  info 
pitch  angle,  heave  acceleration  into  draft  and  the  two 
cross  gains.  Using  typical  gains  calculated  in  the  time 
domain  program  developed  in  Chapt  III  the  roots  of  these 
polynomials  are  located  and  it  is  shown  that  pitch 
dynamics  may  be  very  closely  approximated  by  a 
second-order  system.  Heave  dynamics  are  also 
simpli tied. 

Using  a simplified  approach  to  the  frequency  response 
of  the  craft  to  non  calm-water  conditions,  the  basics 
are  developed  for  later  sea-state  model  worit.  A dynamic 
system  modeling  program  available  at  tie  Naval 
Postgraduate  School  is  used  to  generate  Bode  plots  for 
the  craft’s  frequency  response  in  Heave  and  Pitch  under 
the  simplifying  assumptions. 

3.  SIGNAL  FLOW  GRAPH  ANALYSIS 

Fig  13  was  generated  using  the  sensitivity 
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coefficients  detined  ir.  Jhapter  III.  Fro*  that  griph  the 
following  feed-forward  »al  feedback  loons  and  piths  were 
ident it ied : 


1.  Feed  Forward,  J ( S r A)  intj  r H E T A . 

PI  = (1/IYY)*  (>/S**2) 

2.  Feed  Forward,  J (2)  m*  a Z. 

P 2 = ( 1 Ad  ASS)  « ( 1/S«*2) 

3.  Feed  Forward,  U ( H ? r A ) into  Z. 

? 3 = ( 1/IYY)  *D2TH*  (1/S*«4) 

4.  Feed  Forward,  U(Z)  into  TH3r\. 

P 4 = ( 1/d  ASS)  *D  T HZ  * ( 1/S  * * 4 ) 

5.  Feedback  Loops. 

L 1 = D r HT H/S  * *2 

L2  = DTHDTH/S 

L3  = ( DZTH ) * (DTHZ)  /S**4 


k 

i- 

i 

f 
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L4  = ( D 2d  B ) * (0dB2)  /i  **  3 


L 


L 5 = DZZ/S**2 

Lo  = dhbhb/s 

I 

L 7 = (DZPH)  * (DSBZ)  * (DrBiB) /S**5 


C.  MASON'S  SAIN  SULE  APPLIED 

i i 

M 

.ldson'3  gain  rule  for  a closed  multi-loop  system  is: 

1 

I = SUM  (?k  *DELTA!c)  /DELIA 

j 

where  , 

I = tra nsmiss ion  gain  from  an  input  to  an  output 
node. 

DELTA  = the  graph  determinant 
DELTA*  = cofactor  of  the  *th  path 

P*  = gain  of  the  direct  path  from  input  to  output 

t 

Using  Mason's  gain  rale  the  graph  determinant  (DELTA) 
is  found  to  be  of  the  form: 

* 

D (s)  = (3  5 *S  **5  ♦ 94*s*«4  «■  3 3 *s  **  3 +•  B2*s*»2  + B1*s  * 

SO) /S*  * 5 

The  forward  direct  path  U (THETA)  to  THETA  is: 

N1(s)  = (1/IYY)  * (1/S**5)  * ( A3*s**3  ♦ A2*S**2  + A1*S  ♦ AO) 

The  forward  path  rJ  (Z)  to  Z was  found  to  be: 
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N 2 i s)  = (1/MASS)  * (1/S**5)  *(A3*s**3  ♦ A2*s  * + 2 + \l*s  ♦ 
AO) 


The  forward  direct  patn  U(THETA)  to  Z is: 

N 3 (S)  = ( 1/IYY)  *DZTH*(S  A0)/S**4 

The  forward  path  for  draft  coupling  into  pitch  (0(E) 
to  THETA  is: 

N 4 = ( 1/MASS)  *DTH2*  (S  * A0)/S**« 


0.  SOLUTION  OF  THE  INDIVIDUAL  POLYNOMIALS 

The  gains  used  below  were  taken  from  a typical  time 
domain  run  as  developed  in  chapter  IV. 

DZ  Z = -2178.0 
DOTH  = -121 . 3 
DZMB  = -5499. 0 
DTHi  = 28.78 
DTHTH  = -22.31 
DTHMB  = 72.1 
DMBZ  = -23.79 
DMBMB  = -73.11 

1 . Solution  of  the  Sc  12k  P^termin  ant 

Inserting  the  nunerical  values  for  the  typical 
sensitivity  coefficients  listed  above  and  using  only  the 
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dominant  terms  where  possible,  the  roots  of  the  graph 
determinant  were  found  to  be: 

n • 

DELTA  = (S**2  + 7 1 . 8 *S  + 2066)*(S  + 0 . 5 086 ) * ( 5 * * 2 ♦ 

0. 72*S  ♦22.74) /S**5 

These  factors  represent  a pair  of  complex  poles 
located  at  -35.71+  j27.34  , a real  pole  located  at 
-0.5085  , a pair  of  complex  conjugate  roots  located  at 
-0.3o0  + j 4.755  , and  a fifth  order  aero  located  at  the 
origin . 

2.  Solution  of  Ml  (S) 

N 1 (S)  = (1/IYY)*(S  + 0.  427)*  (S**2  + 72.68*S  + 2l4b)/S**5 

This  polynomial  factors  into  a real  zero  located 
at  -0.427  , a pair  of  complex  conjugate  zeros  located  at 
-3o . 34  + j 28.74  , and  a fifth  order  pole  located  at  the 
origin . 

3 . Solution  of  IU$L 

N 2 ( S)  = ( 1/MASS)  * (5**2  + 0.7^21+S  + 22.3o)*(S  ♦ 

72.38)  /S**5 

which  factors  into  a pair  of  complex  conjuaate 
zeros  located  at  -0.335o  + j 4.712  , a real  zero  at 

-72.38  , and  a fifth  order  pole  located  at  the  origin. 
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4 . Solution  of  n 3 (S| 


N3(S)  = ( 1/IYY)  *DZTH*(S  4 73  . 1 1)  /S**4 

M hich  represents  i real  zero  at  -73.11,  and  a 
fourth  order  pole  at  the  origin. 

5 . Solution  of  Nu  (S) 

N 4 (5)  = (1/MASS)  *DrH^^(S  4 73.11) /S*  *4 

Which  has  the  sane  pole  and  zero  locations  as 

N3 (5)  . 


S.  SIMPLIFICATION  OF  TRANSFER  FUNCTIONS 


ii 


’•  111*1 


T1  (S)  = THETA  (5)  /U  (THETA,  S3 

Therefore  from  Mason's  Sain  Rule: 

1 1 (S)  = Ml  (5)  / DELTA 


T1(S)  = (1/IYY)  * ( S 4 35  . 34  * 428.  7)  * (S  <■  0.427)/ 

(S  35.1  ♦ j 27 . 84)  * (S  «•  0.506)  * (S  ♦ 0.3b  ♦ j 4.  755) 
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To  simplify  let  the  complex  conjugate  zeros  at 
-Jtj.lu  + j2d.7  cancel  the  poles  at  -33.9  + j27.8u  , and 
the  real  zero  at  -3.427  cancell  the  pole  at  -0.508. 
The  result  is  the  simplified  second  order  system: 

T 1 ( S)  = ( 1/IYY)  / (5**2  *■  0 . 72*s  ♦ 22.74) 

Inspection  of  the  time  domain  pitch  response  with 
respect  to  natural  fregueny  and  damping  confirms  this 
basic  underiamped  second  order  response. 


r 2 (S)  = 2 (S)/U  (2,3) 


T 2 (i)  * (1/3ASS)MS  + 0.3356  + j4.712)*(S  ♦ 72.  33)/ 
(S**2  *-35.  1 + j 27.34)*(s  + 0.5035)*((S  +0.33  + 
j4.755) 


Canceling  the  pole  at  -0.3b  + j 4.755  with  the 
zero  at  -0.3S5o  + ju.712  results  in  the  simplified 

z *5  s po 

r2{5)  = { 1/MA3S)  * (S  ♦ 7 2 . 3 3)  / ( 3 + 0 . 508)  * (3**  2 + 71.9*S 

* 2066) 


Further  simplif  ication  can  oe  made  if  (S  + 
72.33)  / tS  ♦ 0.503)  is  approximate!  by  72.33/(S  + 0.503) 
(which  is  valid  for  longer  time  period  considerat ions. ) 
As  a result  T2(S)  becomes: 
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T2(S)  = (7:.  33/HASS)/($0.  508)  *(3«*2  •-  7).8*S  * 2066) 

wnich  has  poles  locate!  at  -0.508,  and  -J5.  4*  j27.33 

This  is  in  quite  close  agreement  with  the  results 
obtained  by  Thaler  and  3erba  14  for  the  same  airflow 
condition  (Kq  * 1.0)  for  their  aeave-only  linearized 

model.  Their  results  were  a real  pole  at  -0.33  and  a 
complex  pair  at  -35.  U 134.4  . This  is  of  the  form: 

T(o)  = K/(S  ♦ 0.38)*(S**2  * 70.2*5  * 24  1 5) 


3 • 11151 

T3  (5)  = 2 ( S)  /U  (THETA  ) 

Here  it  is  noted  that  there  are  no  pole  zero  pairs 
that  conveniently  cancel  so  that  toe  transfer  function 
becomes: 

T3(S)  = { l/IYY)  *D2TH*5*  (5  * 72.33)/ 

(5  ♦ 35.1  £ j 27.84) « (St  0 . 5085 ) * (S  «■  0.3o  ♦ j a . 7 75) 

Making  tae  simplification  used  m the  previous 
section,  the  result  is; 


T3  (5)  = (72.33/IYY)  *32  2 3*5/(5**2  7). 8*5  «•  2055)* 

(S  ♦ 0 . 508)  * (S**  2 t 0.72*5  t 22.74) 

Which  can  be  furtaer  simplified  oy  loosing  at  the 
long  teca  response  vice  initial  few  milli-seconds  to 
(S/ (S  ♦ 0.508))  = 2 so  that: 
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13(3)  = ( 1 4 2/1 Y It ) * DZTH  / (5**2  ♦ 7l.3*s  ♦ 20oS)*(-**2  ♦ 

0 . 7 2*s  ♦ 22.  74) 


a.  ilisi 


* 

U 


Since  the  transfer  fundtions  are  the  sane  except 
for  the  gains  the  solutions  are  the  sane  except  for  the 
gains : 

14(3)  = ( 1/MASS)  *DIHZ  * 5*(S  *■  72. 33) /DELTA 


R 


where  the  same  simplifying  assumptions  may  be  applied. 


?.  SIMPLIFIED  30DE  PLOT  ANALYSIS 


1.  Simplifying  Assumption 
(see  fig.  12) 

The  assumption  is  made  that  the  craft  linearized 
frequency  response  transfer  functions  can  be  separated 
into  two  transfer  functions. 

The  first  block  G1(r4,D,...)  expresses  the 
disturbing  moments  and  accelerations  generated  by  the 
craft  interacting  with  the  waves.  This  function  was 
developed  by  Thaler  and  Gerba  D 4j  for  the  heave-only 
model.  In  that  simpler  model  the  function  31  was  a 
complex  function  of  encounter  frequency,  ingle  of 
incidence,  craft  length  compared  to  encounter  wavelength 
etc..  This  portion  of  the  overall  frequency  response  of 
this  model  is  left  for  future  work  due  to  time 


»i 


constraints. 


The  second  block  32  (We)  represents  the  heave-pitcn 
dynamics  of  the  model  is  developed  in  chapters  II,  III, 
and  IV,  which  respond  to  the  disturbing  moments  and 
accelerations  from  block  31. 

The  Transfer  functions  T1  and  T2  (unsiaplif ied) 
were  simulated  using  the  I BM-3  bO  simulation  lan gu age  DSL 
with  S = 1W  with  the  following  results. 

2.  Bode  Pipt_3  for  Pitch  Dv  na  arc  s 

(see  figs.  15  and  16  ) 

Figure  15  represents  the  normalized  (numerator 
gain  constants  set  to  1|  gain  in  d3's  versus  Omega. 
Note  that  the  abscissa  scale  is  Omega  in  powers  of  10 
since  this  is  a convenient  means  to  achieve  a 
logarithmic  scale  on  the  *-Y  plotter  output.  As  can  be 
seen  the  response  strongly  supports  tne  simplification 
to  second-oruer  system.  These  plots  represent  the  full 
transfer  function  T1  ( 7 W » vice  the  simplified  system. 
Figure  16  further  supports  tne  simplif icati on  as  the 
pnase  plot  starts  at  zeco  degrees  and  shows  only  1 very 
slight  positive  increase  prior  to  the  rapid  change 
toward  -130  characterist ic  of  an  underdamped  system. 
Also  note  that  above  the  damped  natural  frequency  the 
nagnitude  plot  falls  at  -12  dB  per  octave. 

J.  Sou  e Plots  for  Heave  Dv  namics 

(see  figs . 17  and  1 3) 


Figure  17  represents  the  normalized  magnitude 
versus  Omega  far  the  heave  dynamics.  The  system  is 
characteristic  of  an  overdamped  third  order  system  with 
small  gain  increases  at  the  craft's  natural  pitch 
frequency,  which  represents  the  pitch  dynamics  feeding 
back  into  heave  dynamics.  However  the  gain  magnitude 
curve  nas  already  fallen  20  d3  below  it's  lower 
frequency  response  gain  oefoce  the  pitch  coupling  begins 
to  manifest  itself.  As  can  be  seen  in  fig.  18  the  phase 
plot  shows  tne  same  third  order  pole  response  with 
slight  second  order  zero  influence  at  the  natural  pitch 
f requancy . 


VI 
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As  shown  in  Ihapters  TV  ini  V the  simplified 
Pitch-Heave  model  developed  herein  is  a very  reasonable 
simulation  of  the  0-3  test  craft  for  tne  given 
operating  conditions.  Heave,  Pressure,  and  C.G. 
Acceleration  ware  in  vecy  close  agreement.  Pitch  on  the 
other  hand  snowed  some  iiffer^nce  in  natural  frequency 
and  damping  factor.  It  is  felt  that  this  difference  is 
primarily  attriDutable  to  the  inclusion  of  a redundant 
planing  force  calculation  in  the  6 DO?  model, 

historically  the  results  of  the  d DDF  computations  14 
for  pitch  response  have  been  very  similar  to  the  results 
for  the  linear  model,  in  both  natural  frequency  and 
damping  charact er istics . This  can  be  attributed  to  the 
fact  that  the  additional  planing  focce  lid  not  exist  in 
tne  original  5 DO?  program. 


.Since  one  of  the  major  concerns  with  the  o DDF  mod>=i 
is  the  excessive  coaputat ional  times  involved  in 
sea-state  simulations  tne  conclusion  may  be  drawn  that 
eventual  substitution  of  the  sidewall  and  seal  force 
linear  model  equation  in  tne  d DDF  program  will 
substantially  reduce  G?J  time,  with  some  resultant  loss 
in  computational  accuracy.  By  way  of  comparison, 
disregarding  compiling  time  and  output  time,  the  linear 
model  required  approximately  9 seconds  of  CP'J  time  for 
the  5 second  simulation  while  the  o DDF  model  required  9 
minutes  of  c?J  for  the  5 second  run.  This  is  a 
computation  time  improvement  on  tne  order  of  60  times. 

From  the  time-domain  plenum  gauge  pressure  data  it  can 
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oe  seea  that  the  dynamic  variation  in  effective  plenum 
roof  area  has  a noticeable  effect  on  the  plenum  pressure 
transient  behavior  as  compared  to  previous  resuLts  14  . 
The  initial  pressure  down  transient  was  faster,  dropped 
to  a lower  value  , and  began  it's  recovery  more  guickly. 
Therefore  the  pLenua  roof  dynamic  variations  strongly 
affects  the  C.G.  acceleration  characteristics  of  the 
craft  and  therefore  it* 3 transient  behavior. 


The  dynamic  center  of  pressure  variation  developed 
herein  proved  to  be  a very  draft  sensitive  phenomenon. 
Although  the  length  of  the  other  lever  arms  as  modeled 
herein  change  by  a similar  amount,  percentage  of 
variation  is  radically  different.  For  the  seals  the 
variation  is  on  the  order  of  inches  in  10  feet,  where 
the  variation  in  lever  arm  for  plenum  pressure  lift  is 
on  the  order  of  inches  in  a single  foot  (or  less.)  This 
implies  that  dynamic  3.P.  variation  is  the  dominant 
factor  in  steady  state  pitch  angle  versus  steady  state 
draft. 


Planing  forces  as  included  in  this  model  are  now  in  a 
form  where  the  basic  factors  as  far  as  craft 

construction  characteristics  are  readily  recognized. 
This  force  provides  part  of  the  answer  to  craft 

pitch-down  characteristics  as  initially  investigated  by 
Reidel  C91  , with  the  remainder  of  the  moment  being 
developed  by  the  C.P.  variation 
planing  force  provides  a clue 
construction,  the  craft’s  pitch 
character istics  may  be  controlled. 


for  pressure.  The 
as  to  how,  during 
amplitude  excursion 


Several  areas  for  future  investigation  have  become 
apparent  in  this  work. 


1.  Verification  of  the  effect 


of  dynamic  center  of 


pressure  variation  and  planing  forces  as  incorporated  in 
the  6 DOF  model  over  the  entire  speed  operating  range. 


2.  Incorporate  the  simple  modeling  for  sidewalls  and 
seals  into  the  6 DOF  modal  and  investigate  the  trade-off 
between  CPU  time  and  computational  accuracy. 

3 . A closer  look,  at  the  seal  forces  used  herein.  It 
is  possible  that  a simple  gain  factor  needs  to  be 
included  in  these  eguations  in  order  to  optimize  the 
steady-state  seal  forces  with  respect  to  actual  load 
data  as  measured  during  experimental  test  runs  of  the 
craft  10  . 

4.  Incorporate  sea  state  forces  and  moments  into  the 
simplified  model  by  simulation  of  the  forces  and  moments 
that  are  generated  by  craft  interaction  with  the  sea. 
Evaluate  the  accuracy  of  the  method  used  by  comparison 
with  data  generated  by  the  simple  model  and  the  6 DOF 
model  in  sea-state  conditions. 

5.  Analytical  investigation  of  the  craft  physical 
characteristics  that  have  predominant  influence  on  the 
natural  pitch  frequency  and  damping  oy  using  the  gain 
terms  from  the  results  of  the  signal  flow  grapn  analysis 
done  in  Chapter  V. 

6.  Use  the  simplified  modal  to  develop  a ride-comfort 
control  system.  lest  the  resulting  design  on  the  6 DOF 
aodel,  and  compare  the  results  ootained  with  those 
suggested  by  3D  HR  C7]  . It  mas  oeen  snown  empirically 
for  some  time  that  the  era  ft  ex  .limited  simple  second 
order  behavior  in  the  pitch  mode.  Mow  that  this  fact 
has  been  analytically  verified  the  implications  for  an 
automatic  control  system  for  pitch  response  are 
powerful.  A simple  control  system  can  be  visualized 
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that  does  not  involve  either  complex  optimal  control 
system  design  work  or  hardware.  \ simple  acceleration 
and  rate  feedback  contral  system  is  implied,  with  the 
possibility  of  adding  a speed  input  to  the  system  to 
compensate  for  the  localized  linearity  of  the  craft 
character i sties. 
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LINEAR  SES  MODEL  USER'S  GUIDE 


A.  General  Comments 


The  program  as  listed  following  Appendix  3 is  intended 
for  useage  via  the  hot  card  reader  at  NPS. 

It  is  intended  to  run  in  Foctraa  G with  Plotting 
option  (PROC  = EOP.TCLGP.) 

The  program  is  structured  on  the  assum-ption  that  the 
initial  and  final  state  variables  are  in  equilibrium 
prior  to  data  entry,  but  will  run  regardless.  The  force 
and  moment  residual  listings  will  indicate  the  value  of 
imoalance,  if  any. 

The  craft  physical  dimensions  are  alterable  by 
changing  the  cards  in  the  craft  dimension  section  of  the 
main  program,  which  are  keyed  to  the  craft  geometry 
nerein . 

3.  Data  Deck  Setup. 

The  /*  and  //GO.sysiN  DD  * cards  are  inserted 
following  the  end  card  in  Suoroutine  VERSA?,  followed  by 
the  data  deck  as  described 

1 .Plot  option  card  (TPLOT)  , format  IT.  IPLOT  = 1 = 

VSRSATEC  X-y  plotter  output. 

TPLOT  = 0 = PST  PLOT  output 
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i 

. 4 


2.  List  option  card  (ILIST),  Format  II 

ILISI  = 1 = Forces  and  moments, with  residuals  and  all 
sensitivity  coefficients  listed. 

ILIST  = 0 = Input  craft  operating  conditions  only 
listed . 

3.  Time  and  Integration  variables. 

TI,  TF,  NPTS,  NSTSP  Format  (2?10.4,2I10) 

Ti  = problem  start  time  (normally  0.0) 

TF  = Problem  stop  time  (seconds) 

NPTS  = number  of  tabular  output  points  desired. 

NSTEP  = number  of  integration  steps. 

note:  it  is  best  to  mafce  NSTSP  an  integer  multiple 

of  NPTS , and  never  less  than  NPTS. 

4.  Planing  coefficient  and  arbitrary  center  of  pressure 
input,  i.  e.  aSAO  ( ) AKP,  PLCOEFF 

Form  at  = 2F 1 0 . 4 

5.  Initial  Conditions  card. 

The  initial  craft  conditio  os  ace  read  in  in  the 
following  order:  weight (Lbs .) , Draft  (in).  Pitch  angle 
(Deg.)),  F53AS  (psf),  and  speed(fps.)  Format  is  5F10.4 

6.  Final  Conditions. 

The  final  conditions  are  read  in  the  following  order 
and  format:  weight  (Lbs),  Draft  iin).  Pitch  angle 
(Deg),  and  P3BAR  (psf)  , Format  4F10.4.  Since  this  is  a 
constant  speed  model,  initial  and  final  speed  are 
assumed  equal. 

7.  Label  Cards. 


Label  cards  for  VE3SAIEC  output  are 


he  last  10  carls. 


There  are  five  graphs  in  this  node.  Draft(in),  C. G. 
Acceleration  (G's)  , Pitch  Angle  (Deg.)/  Pitch 
Acceleration  (Deg/SEC* *2)  , and  P33A8  (Each  graph 
reguires  two  label  cards.  Format  (bA8) , which  print  on 
two  lines  below  the  graph.  The  graphs  are  output  in  the 
order  listed  above. 


1 

; 
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APPENDIX  3 


KEEL  CROSS  SECTION 


FIGURE  3A 


VERTICAL  GROSS  SECTION 
FORWARD 

FIGURE  3B 


VERTICAL  CROSS  SECTION 
AFT 


FIGURE  3C 


Figure  3 - ASSUMED  GEOflETS'f  (CONT* 
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X-SCPLE= 1 . OOE+GG  UNITS  INCH. 

T -SCALE -5 . GOE +GG  UNITS  INCH. 

plot  is  plenum  gouge  - ~ 

COB  SES  LINEAR  2 OQ'r  V1G 

figure  5 - LINEAR  PLENUM  PRESSURE  TRANSIENT 


X -5CRLE  = 1 . 0 0 E G C UNITS 
T-SCRLE=2. 00E-G2  UNITS 

PLGT  IS  CG  RCG 
CflS  SES  LiNh^R 


NCH. 

NCM. 


2 D C 


m g : 


Figure  7 - LINEAR  C.J.  AC  CEL  £R  AT  ION 


1 


03 


X-SCfiLE=l . 00E-00  UNITS  INCH. 
T-SCflLE=l . 00E-01  UNITS  INCH. 

PLOT  IS  PITCH  HNG 
CRB  SES  LINEAR  2 

Figure  3 - LINEAR  PITCH  TR 


IlJ  O'" 


AD-A070  368 


UNCLASSIFIED 


NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CA  F/G  13/10 

DEVELOPMENT  AND  TIME  DOMAIN  VALIDATION  OF  A LOV-ORDERf  CONSTANT— ETC (U) 
MAR  79  L V BARNES 


000  001  002  003  00U  005 


X-SC ALE= 1 . OOE+OO  UNITS  INCH. 
T-SCRLE=2. OOE+OO  UNITS  INCH. 

30  KNOT  RUN 

PLOT  IS  Z DISPLACEMENT 

Figure  10  - 6 DOF  DB AFT  TRANSIENT 
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X-SCRLE= 1 . 00E+00  UNITS  INCH. 
Y-SCRLE=2. 00E-02  UNITS  INCH. 

30  KNOT  RUN 

PLOT  IS  C, Go  ACCELERATION 

Figure  11  - 6 DOF  C.G.  ACCELERATION 
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X-SCALE= l . OOE+OO  UNITS  INCH. 

' -SCALE* l . 006 - 0 1 UNITS  INCH. 

30  KNOT  RUN 
PLOT  15  PITCH  ANGLE 

\2  - b DOF  PITCH  TRANSIENT 


Pi  5 jr* 


XSC  ALE  = 0 . SO  UNITS/INCH  RUN  NO.  1 

y sc  ale  = io.og  units/inch  plot  no.  : 


Figure  15  - LINEAR  PITCH  30DE  PLOT  t N \ 3 > 


00*0*1  00*0  OO'Ob-  00*00-  OO'Oei-  00*031-  00*00?- 


ou 


X5CRlE=  0.80  UNITS/INCH  RON  NO.  1 

T SCRlt  = 40.00  UNITS/INCH  PlOT  NO.  0 

i 

io  - k.nun  tfiiLn  ojub  plUT  *■  fc»rl  A 3 E ' 

* * 
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1XSC ALE=  0.80  UNITS/INCH  RUN  NO.  1 | 

T5CflLE=  3.00  UNITS/INCH  PLOT  NO.  1 

*igur«i  17  - LTME1H  HEAVE  30DE  PLOT  < H A i > 
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i XSC ALE=  0.80  UNITS/INCH  RUN  NO.  1 

T SCALE  = 40.00  UNITS/INCH  PLOT  NO.  2 

Figure  18  - LINEAR  HEAVE  BODE  PLOT  <PHASE> 
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CALL  VERSAP  (N  PTS  , TI  M E , PLOT  V) 

GO  TO  3 5 

44  CALL  PLOTP (TIME, PLOTV, NPfS, 0)  SES03260 
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